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Abstract. We consider a stochastic differential equation involving standard 
and fractional Brownian motion with unknown drift parameter to be esti- 
mated. We investigate the standard maximum likelihood estimate of the drift 
parameter, two non-standard estimates and three estimates for the sequential 
estimation. Model strong consistency and some other properties are proved. 
The linear model and Ornstein-Uhlenbeck model are studied in detail. As 
an auxiliary result, an asymptotic behavior of the fractional derivative of the 
fractional Brownian motion is established. 

[2010]60G22; 60J65; 60H10; 62F05 Fractional Brownian motion, Brownian mo- 
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1. Introduction 

Modern mathematical statistics tends to shift away from the standard statistical 
schemes based on independent random variables; besides, these days many statis- 
tical models are based on continuous time. Therefore, the corresponding statistical 
problems (e.g., parameter estimation) can be handled by methods of the theory 
of stochastic processes in addition to the standard statistical methods. Statistics 
for stochastic processes is well-developed for diffusion processes and even for semi- 
martingales (see, for instance, |LipSh| ) but is still developing for the processes with 
long-range dependence. The latter is an integral part of stochastic processes, fea- 
turing a wide spectrum of applications applications in economics, physics, finance 
and other fields. The present paper is devoted to the parameter estimation in such 
models involving fractional Brownian motion (fBm) with Hurst parameter H > ^ 
which is a well-known long-memory process. The paper also studies a mixed model 
based on both standard and fractional Brownian motion which turns out to be 
more flexible. One of the reasons to consider such model comes from the modern 
mathematical finance where it it has become very popular to assume that the un- 
derlying random noise consists of two parts: the fundamental part, describing the 
economical background for the stock price, and the trading part, related to the 
randomness inherent to the stock market. In our case the fundamental part of the 
noise has a long memory while the trading part is a white noise. 
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Statistical aspects of models involving fractional Brownian motion were studied 
in many sources. One of the important problems in particular is the drift parame- 
ter estimation. In this regard, let us mention papers 'HuNu' and [ KlLeBrj . where 
the fractional Ornstein-Uhlenbeck process with unknown drift parameter originally 
was studied, books [BishOS] . [MishOSj and |Prara| and the references therein, and 
papers [BTTj . [XZXj . [XZZj . and [HuXZ I , where the estimate was constructed via 
discrete observations. We shall also use the results for sequential estimates for 
semimartingales from |MN88] . In the present paper we consider stochastic differen- 
tial equations involving fractional Brownian motion along with equations involving 
both standard and fractional Brownian motion. We derive the standard maximum 
likelihood estimate and propose non-standard estimates for the unknown drift pa- 
rameter. Several non-standard estimates for the drift parameter were proposed in 
[HuNu] for the fractional Ornstein-Uhlenbeck process. We go a step ahead and 
propose non-standard estimates for the drift parameter in a general stochastic dif- 
ferential equation involving fBm. For the models involving only fractional Brownian 
motion, we compare properties of the estimates. In the mixed models the standard 
maximum likelihood estimate does not exist but the non-standard estimate works. 
To formulate the conditions for strong consistency of the non-standard estimates, 
we need to investigate the asymptotic behavior of the fractional derivative of the 
fractional Brownian motion using the general growth results for Gaussian processes. 

The paper is organized as follows. In Section 2 we introduce the models and 
the estimates: the maximum likelihood estimate, two non-standard estimates and 
three sequential estimates. Asymptotic growth of the fractional derivative of fBm 
is established in Section 4. Section 5 contains the main results concerning the 
strong consistency of all estimates and some additional properties of sequential 
estimates. The linear model and Ornstein-Uhlenbeck model are studied in detail. 
We generalize the result of strong consistency of the drift parameter estimate in the 
Ornstein-Uhlenbeck model from [KlLeBrj to the model with variable coefficients. 

2. Model description and preliminaries 

2.1. Model description. Let {il, J^, P) be a complete probability space with 
filtration J- = {J-j, t € M"*"} satisfying the standard assumptions. It is assumed that 
all processes under consideration are adapted to filtration 

Definition 1. Fractional Brownian motion (fBm) with Hurst index H e (0,1) is 
a Gaussian process = {B^,t G R+} on (i7, J^, P) featuring the properties 

(a) B^=0; 

(b) EB^ = 0,t€R+; 

(c) EBl^B^ = i(i2^^ + _ |t _ s\^"),s, t e M+. 

We consider the continuous modification of B^ whose existence is guaranteed 
by the classical Kolmogorov theorem. 

To describe the statistical model, we need to introduce the pathwise integrals 
w.r.t. fBm. Consider two non-random functions / and g defined on some interval 
[a, b] C M"*". Suppose also that the the following limits exist: f{u+) := lim^j^o f{u + 
5) and g{u—) := lim^^o g(u — (5), a < u < b. Let 



fa+{x) := (/(a;) - f{a+))l(^a,b){x), gb-{x) := {g{b-) - g{x))l(^a,b){x). 
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Suppose that fa+ G /^^_(Lp[a, 6])), gt- G llz°'{Lq[a,b])) for some p > 1, q > 
l,l/p+l/q < 1, 0<Q!<1. (For the standard notation and statements concerning 
fractional analysis, see |SMK| ). Introduce the fractional derivatives 

1 / /a+(s) , r fa+{s) ~ fa+{u) 



iT^oc f \l \ ^ ( Ja+(S) / Ja+(S ~ Ja+(U \ , , 



It is known that 'D^^fa+ G Lp[a,b], Vlz'^gb- G Lg[a, 6]. 

Definition 2. ( |Zah98| . |Zah99| ) Under above assumptions, the generalized (frac- 
tional) Lebesgue-Stieltjes integral f{x)dg{x) is defined as 

[\K+fa+){x)iVlz''gb-){x)dx + f{a+){g{b-) - g{a+)), 
and for ap < 1 it can be simplified to 

'f{x)dg{x) :=e-" f [VZ^f){x){Vlz'' gb-){x)dx. 



As follows from [SMKj . for any 1 — H < a < 1 there exist fractional deriva- 
tives Vlz°^B^_ and Vlz°'Bl^_ G Loo [a, 6] for any < a < 6. Therefore, for 
/ G /"_|_(Li[a, 6]) we can define the integral w.r.t. fBm in the following way. 

Definition 3. f [NuaR| . |Zah98| . |Zah99j ) The integral with respect to fBm is de- 
fined as 

(1) f fdB" := f\v:^f){x){Vlz-Bt){x)dx. 

J a J a 

An evident estimate follows immediately from ([1]): 

(2) f fdB" < sup \{Vlz^'B^_){x)\ [' \iV:^f )ix)\dx. 

Ja a<x<b Ja 

Let US take a Wiener process W = {Wt,t G M+} on probability space {ft, J-, T , P), 
possibly correlated with B^ . Assume that H > ^ and consider a one-dimensional 
mixed stochastic differential equation involving both the Wiener process and the 
fractional Brownian motion 



(3) Xt = XQ + 



' [ ais,Xs)ds+ [ b{s,Xs)dB^+ [ c{s,Xs)dWs, i G R+, 
Jo Jo Jo 



where G M is the initial value, 9 is the unknown parameter to be estimated, the 
first integral in the right-hand side of (jS]) is the Lebesgue-Stieltjes integral, the sec- 
ond integral is the generalized Lebesgue-Stieltjes integral introduced in Definition 
[HI and the third one is the Ito integral. From now on, we shall assume that the 
coefficients of equation ([3]) satisfy the following assumptions on any interval [0,T]: 
(^i) Linear growth of a and b: for any s G [0,T] and any a; G K 

\ais,x)\ + \b{s,x)\<K{l + \x\). 
{A2) Lipschitz continuity of a,c in space: for any t G [0,r] and G M 
\a{t, x) - a{t, y)\ + |c(t, x) - c(i, y)| < K\x - y\. 
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{A3) Holder continuity in time: function h{t,x) is differentiable in x and there 
exists /3 G (1 — i/, 1) such that for any s, t € [0, T] and any a; g R 

|a(s, x)-a{t, x)\ + \h{s, x)~b{t, x)\ + \c{s, x)-c{t, x)\ + \d^b[s, x)-dxb{t, x)\ < K\s-tf . 
(Ai) Lipschitz continuity of dxb in space: for any t € [0,T] and any a;, y G M 

\d,b{t,x)~d.M,y)\<K\x-y\. 
{Aq) Boundedness of c and dxb: for any s G [0,r] and x G M 

\c{s,x)\ + \dxb{.s,x)\ <K. 

Here K is a, constant independent of a;, y, s and t. For an arbitrary interval 
[0, T], a > and k = ^ A /3 define the following norm: 



ll/lloo,a,[o,T] = sup {\fis)\+ |/(s)-/(z)|(s-z)-^-"dz 

s£[0,T] \ Jo 

It was proved in |MiSh) that under assumptions (Ai) — {Ac,) there exists solution 
X — {Xt,J^t,t G [Oi^]} for equation ^ on any interval [0,T] which satisfies 

(4) ll^lloo,a,[0,T] < 00 a.S. 

for any a € {I — H, k). This solution is unique in the class of processes satisfying 
(U for some a > 1 ~ H . 

Remark 1. In case when components W and are independent, assumptions for 
the coefficients can be relaxed, as it has been shown in jGuNuj . More specifically, 
coefficient c can be of linear growth, and dxb can be Holder continuous up to some 
order less than 1. 

2.2. Construction of drift parameter estimates: the standard maximum 
likelihood estimate. To start with, consider the case c(t, a;) = which was stud- 
ied, for instance, in |KlLeBr) and |Mish08] . Recall some facts from the theory of 
drift parameter estimation in this case. Consider the equation 

(5) Xt^xo + 9 [ ais,Xs)ds+ [ b{s,Xs)dBf, i G M. 

Let assumptions (^i) and {A^) with c = hold on any interval [0,T], together 
with the following assumptions: 

{A'2) Lipschitz continuity ofa,b in space: for any t G [0,r] and a;,y G M 

\a{t, x) ~ a{t, y)\ + a:) - b{t, y)\ < K\x ~ y\, 

{A'^) Holder continuity of dxb{t,x) in space: there exists such p G (3/2 — iJ, 1) 
that for any t G [0, T] and x,y 

\dxbit,x)-dxb{t,y)\ < Dlx-yl", 

Then, according to [NuaRj . solution for equation ([5]) exists on any interval [0,T] 
and is unique in the class of processes satisfying Q for some a > 1 ~ H. 
In addition, suppose that the following assumption holds: 

(Bi) b{t,Xt) ^ 0,t G [0,r] and is a.s. Lebesgue integrable on [0,T] for 

any T > 0. 
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Denote il^{t,x) = , ip{t) := il;{t,Xt). Also, let the kernel 

with ch — ( 2HT(^^H)^(H+ j- ) ) ' ^"'^'^ introduce the integral 

(6) Jt= / lH{t,s)ip{s)ds^CH f {t~ s)^-"s^-"ip{s)ds. 

Jo Jo 

Finally, let = Init, s)dB^ be Gaussian martingale with square bracket 

{M)f = t2-2H (Molchan martingale, see INVV99I). 
Consider two processes: 



Yt^ [ b-\s,Xs)dXs ^0 [ ip{s)ds + Bf 
JO Jo 



and 



Zt= f lH{t, s)dYs = 9Jt + Mf . 
Jo 

Note that we can rewrite process Z as 

t 

lH{t,s)b-\s,Xs)dX,, 

so Z is a functional of the observable process X. The following smoothness condition 
for the function ip (Lemma 6.3.2 |Mish08| ) ensures the semimartingale property of 
Z. 

Lemma 1. Let ip{t,x) G Ci(M+) x C2(M). Then for any t > 
(7) 

/■* 

,1-2H 



J'{t) = (2 - 2H)Ch^{0, xo)t^-^" + / lH{t, s) {^'t{s, Xs) + e^'^is, Xs)a{s, X,)) ds 

Jo 

^{h-]^ch s-^-"{t-s)^-" (ip[{u,X^) + ei>'^{u,X^)a{u,Xu))duds 

+ {2 - 2H)cHt^-^" f s^"-"" f ui-"{s-uf2-H^'^(u,X^)b(u,Xu)dB^ds 
Jo Jo 

+CHt-^ I ui-"{t~u)i-"^'Ju,Xu)b{u,Xu)dB^, 



where Ch = B{^ ~ H,^ — H)ch = 2 gr ( g ) rx 3 - 2 h ) ) ' '^"''^ "^^^ '^^ involved 
integrals exist a.s. 

Remark 2. Suppose that e Ci(M+) x C^{R) and limit ^(0) = lim3^o<r(s) 

exists a.s., where <^(s) = S2^^(p(s). In this case J{t) can be presented as 

J{t)^CH / (t-g)^-^^(,s)d5^ f ' m + CH / ^ 3 ^ ^ <,'{s)ds, 
JO 2 ~ J 2 ^ ^ 
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and J'{t) from ^ can be simplified to 

+0i^Us,Xs)a{s,Xs))ds+ f lH{t,s)^P',{s,Xs)b{s,X,)dB^. 



Same way as Z , processes J and J' are functionals of X. It is more convenient 
to consider process = (2 - 2H)-^ J' {t)t^"-^ , so that 

Zt = {2~ 2H)9 f xis)s^-'^"ds + Mf ^0 f x(s)d(Af + Mf . 
Jo Jo 

Suppose that the foUowing conditions hold: 
(B2) EIt ■■= E xld{M"), < 00 for any T > 0, 
(S3) /oo := J^xld{M"), = 00 a.s. 
Then we can consider the maximum likelihood estimate 

,(1) _ loXsdZ, ^ ^ loXsdM^ 



Condition {B2) ensures that process XsdM^,t > is a square integrable mar- 
tingale, and condition (i?3) alongside with the law of large numbers for martingales 

ensure that Ji? ^° — ^— — > a.s. as T ^ 00. Summarizing, we arrive at the 

Jo Xi<l(MH)^ 

following result ( |Mish08) ) . 

Proposition 1. Letipit.x) G C^(K+) x C^(R) and assumptions {Ai), (A3), {A'2), 
{A\) and (i?i)-(i?3) hold. Then estimate O^rp is strongly consistent as T ^ 00. 

2.3. Construction of drift parameter estimates: two non-standard esti- 
mates. In case when c = 0, it is possible to construct another estimate for pa- 
rameter 9, preserving the structure of the standard maximum likelihood estimate. 
Similar approach was applied in |HuNu| to the fractional Ornstcin-Uhlenbeck pro- 
cess with constant coefhcients. We shall use process Y to define the estimate as 



(8) 



,(2) _ lo ^sdYs ^ ipsdB^ 

rT o , ^ r-T 



lo "fids /o iflds 

Let us return to general equation ^ with non-zero c and construct the estimate 
of parameter 9. Suppose that the following assumption holds: 

(Ci) c{t,Xt) 7^ 0,i e [0,r], is a.s. Lebesgue integrable on [0,T] for any 

T > and there exists generalized Lebesgue-Stieltjes integral /J" c(t'xt) dBt^ ■ 

Define functions ii)i{t, x) = and '02 (t, x) = ^[j^, processes ipi{t) = ipi{t, Xt),i = 

1, 2 and process 

Yt^ f b-\s,X,)dXs^9 f ^i{s)ds+ f <f2{s)dB^ + Wt. 
Jo Jo Jo 

Evidently, y is a functional of X and is observable. Assume additionally that the 

generalized Lebesgue-Stieltjes integral ipi{t)if2{t)dB^ exists and 

(C2) for any T > ipl{s)ds < 00. 
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Denote '&{s) — ipi{s)ip2{s). We can consider the following estimate of parameter 

9: 

(9) e(3) /o'yi(^)^y^ _g , lomdB^ ^ Iofiis)dWs ^ 

vl{s)ds ipl{s)ds J^^ ipl{s)ds 

(3) 

Estimate Orp preserves the traditional form of maximum likelihood estimates for 
diffusion models. The right-hand side of © provides a stochastic representation of 

(3) 

9>p . We shall use it to investigate the strong consistency of this estimate. 

2.4. Construction of drift parameter estimates: sequential estimates. Re- 
turn to model ([5]) and suppose that conditions (Bi) — (B^) hold. For any h > 
consider the stopping time 

' x'd(M). = h}. 



T{h) = inf{i > : / 



Under conditions (Bi) — (^2) we have r(/i) < 00 a.s. and /J^^'*'' x'id{M)s — h. The 
sequential maximum likelihood estimate has a form 



(10) 



h + h ■ 

(2) „„A /i(3) 



Sequential versions of estimates d^p' and 9^' have a form 



and 

/)(3) 

where 



5(2) _/) , io_^fl^ 
" + h 

C^mdB^ , C\^{s)dWs 



t 

2/ 



v{h)^M{t>0: / ipi{s)ds = h}. 
Jo 

To provide an exhaustive study of the introduced estimates, we will need a 
number of auxiliary facts about Gaussian processes. These facts are presented in 
the next section. Technical proofs may be found in Appendix. 

3. Auxiliary results for Gaussian processes related to the 

FRACTIONAL BROWNIAN MOTION. 

We start with the exponential maximal bound for a Gaussian process defined on 
an abstract pscudometric space, expressed in terms of the metric capacity of this 
space. This result is a particular case of the general theorem proved in |BulKoz) . 
p. 100. 

Lemma 2. Let T be a non-empty set, X = {X{t), t G T} be centered Gaussian 
process. Suppose that the pscudometric space (T, p) with pscudometric 

p{t,s)^{E{X{t)-X{s)fY 

is separable and process X is separable on this space. Also, let the following condi- 
tions hold: 

a:^su-p{E\X{t)\^y < 00, 
teT 
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and 

(log Afj'(w))2(iu < oo, 

where Nt{u) is the number of elements in the minimal u-covering of space (T, /o). 
Then for any A > and any 6 G (0, 1) the following inequality holds: 

^exp|Asup|X(<)|l < 2Q{X,9), 
I teT J 

where 

{\2 2 2A Z"^" 1 1 

Consider set T = {t = (fi, ^2) G : < t2 < ii} supplied with the distance 
m(t,s) = |fi-si|V|t2-S2|. 

Assume random process X = {X(t),t E T} satisfies the following conditions. 

{Di) Process X is a centered Gaussian process on T, separable on metric space 
(T,m). 

{D2) There exist /3 > 0, 7 > and a constant C(/3, 7) independent of X, t and s 
such that for any t , s G T 

(11) {E(X(t) - X{s)r) ^ < C{P, 7) (ti V s^f [m (t, s))^ . 

[D^) There exist (5 > and a constant C(5) independent of X and t such that 
for any teT 

(12) {E{X{t)fY <C{5)t{. 

Let us introduce the following notations. Let A{t) > l,t > be an increasing 
function such that A{t) — >■ 00, t — >■ 00. Consider an increasing sequence 60 — 0, 
be, < be+i,l > 1 and suppose that be ^ 00, i ^ 00. For 5e = A{be) and k > we 
denote 

S{6) = £ bi^,6^\ ^, = ^fi + ^-^\ B, = C{6)S{5), 
e=o V 7 7/ 

Ci = C2K~--S{6 + Ki) and C2 = ^^{0(6)^-^ {C{l3,j))^ . 

27 

Now we shall present the auxiliary exponential maximal bound for a Gaussian 
process defined on (T,m). 

Theorem 1. Let {X{t), t e T} be a random process satisfying assumptions (I?i) — 
{D3). LetO<a<b, set Ta,b = {t= {tiM) € T : a < ti < b, < t2 < ti}. Then 
for any < <1, A>0 and < k < 1 A 27 the following inequality holds: 

Eexplx sup \X{t) \ \ < 2Q{X,e), 

[ teTa,b J 

where 

Q(X9)= exp / >^'(b'C(^)r ^ ^ 1 

y(A,y) exp I 2(1-0)2 +i_e'' e^^h]- 
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Proof. It follows from and (HI]) that 

(13) — — fj^\Vf4-\\2\i 



and 
(14) 



d:= sup {E\X{t)\')'' < C{5y 



sup {E{X{t) - X{s)f) ^ < a{h) C(^, -i)hl^h'^. 

m(t,s)</t,t,s GT„_b 



In turn, it follows from that 
h — a 



(15) Nt:^,{v)< 



2a(-i)(w) 



1 



2a(-i)(v) 



1 < 



2v-i 



Define J (6ld) := /p"'* (log ArT„,,(w)) ' du. It follows from ^ that 



(16) 7(6*, 
For any < k < 1, 



i-ed 




/ V2 


log 1^ 


Jo 





(C(/3,7))" ^ 



1 a;'' 
log(l + x) = - log(l + x)'" < — 



Now, let K e (0, 1 A 27). Then it follows from ^ and (HH) that 



J (6lrf) < 

K2 Jo 



V2 r {{c{Pn))^h'+^y: 



■dv 



{2v ) 2 

V2 Ac(/3,7))Hi+f V 



K2 



Separability of X on (T,m) and relation ([T4|) ensure separability of X on (T,p) 

with p(t,s) = (^E{X{t) — X(s))^) ^. Hence the statement of the theorem follows 
from Lemma [5) □ 

Now we are ready to state the general result concerning the asymptotic maximal 
growth of a Gaussian process defined on (T, m). 

Theorem 2. Let X = {X{t),t E T} satisfy assumptions (Di) — (-D3). Suppose 
that function A{t) is chosen in such a way that series S{S) converges. In case when 



1 + 



13 5 



> 0, assume additionally that there exists such < k < 1 that series 



S{5 + Ki) converges with ni = ^ (^1 + ^ — . 

Then there exists such random variable ^ > that on any w £ il and for any 

te T 

< A{t,)t 

and ^ satisfies the following assumption: 
(D4) for any e > (2Ci + 1)^ 

(^s - {2Ci + 1)]^ ^ 



> e} < 2 exp < 



2Bf 



> . 
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Here the value of n < 2^ is chosen to ensure the convergence of series S{6 + ki) in 
case when 1 + - — - > 0, and we set k = i A 7 in case when 1 + - — - < 0. 

Proof. It is easy to check that 

(17) I:=Ee^p\xsnp^-^\<Ee^p\xf2{6i)-' sup |X(t)||. 

I teT ^(iij J [ ^ tie(b«,fc<.+i) J 

Let i > 0, ri > 1 be such integers that J^'eLo ~ ^- Then it fohows from (|T7| . 
Theorem [1] and Holder inequahty that for any 9 G (0, 1) and < k < 1 A 27 

1 

/<nM^exphy sup |X(t)| <n(2Q£(A,^?))^ =2l](Q,(A,(?))^, 

i=Q \ \ tie(b.,6f+i) J / £=0 

where 

C^£(A, 6*) - cxp <^ — 2-— + 7; STi-o 



2<52 (1-0)2 ' {\-B)h"' 
Therefore, if we take such value of k < 27 that series S{6 + ki) converges in case 
when 1 + - — - > and set k = i A 7 in case when 1 + - — - < 0, we obtain 

(18) / < 2exp I ^:(^(^))^ V !#! + 2AC..~^5(. + .1) 1 

^ ^ - 2(1-0)2 J2 (1_0)0#F / 

Now we can substitute = S'((5)&^*(5£ into p^ : 

' A2(S'((5)C((5))2 2AC2K-^5((5 + Ki)l 



2(1-0)2 



/ < 2 exp 

Therefore, 

(19) ^exp/Asupi^^l < 2cxp(^i?2 + 2AC'L 



teT ^K^i) 



A{ti) i -"^^^"12 



where 



B^'Jm^ and c'^^^'^"^^(^ + '^^) 



1-0 (1-0)0- 
It follows immediately from p9|) that for any A > 0, e > 



(20) P (sup ^ > A < exp{-Ae}i?exp (asup ^\ < 

f A2 

< 2exp<^ —B'^ + 2\C-\e 

If we minimize the right-hand side of (j20p w.r.t. A then we obtain that for any 
e>2C 



2 exp ■ 



(£(l-0)-20-^Ci)2 
2Sf 
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Finally, we can insert 6 = e ^t+k into (f21|) and derive that for e > (2Ci + 1)^-1+" 



teT 



Denote ^ := sup ^^rt^- Then ^ satisfies assumption (I?4), and on any a; € f2 

tGT ^ ^' 

X{t) < 

which concludes the proof. □ 
Theorem 3. Let < H < 1,1 - H < a < 1, T = {t = {tiM)-, ^<h< h}, 

r>H r>H ptl F>H tdH 

Then for any p > 1 there exists random variable ^ — ^(p) such that for any t G T 

i^(<)i<((if+"-^(iog(ti)r)vi)e(p), 

where ^{p) satisfies assumption (Di) with some constants Bi and Ci. 

The proof of Theorem [3] is of a technical nature and therefore it is placed in 
Appendix. 

4. Main results 

4.1. General results on strong consistency. In this section we shall establish 
conditions for strong consistency of Oj, and Oj, . 

Theorem 4. Let assumptions [Ai], (A3), {A'2), {A'^) (Bi) and {B2) hold and let 
function ip satisfy the following assumption: 

(-B4) There exists such a > 1 — H and p > I that 

T^+'^-HiogTr lo m+^){s)\ds ^ 

(22) 7p — !• a.s. as J — ^ 00. 

Iq "fids 

12) 

Then estimate is correctly defined and strongly consistent as T —>■ 00. 
Proof. We must prove that ^°^t'^^1 >■ a.s. as T — > 00. According to ([2]), 

r VsdB^ < sup \{V]r^B?_){t)\ f \[V^^^){s)\ds. 

Jo Q<t<T Jo 

Furthermore, according to Theorem[3l for any p > 1 there exists a random variable 
£, = £,{p) independent of T such that for any T > 

sup \{V},-_"B?_){t)\<ap)T"+"~\logTr, 

0<t<T 

which concludes the proof. □ 
Relation (1221 ensures convergence -/o Jf^^-^^ ^ g a.s. in the general case. In a 

Jo nd-'^ 

particular case when function (p is non-random and integral (psdB^ is a Wiener 
integral w.r.t. the fractional Brownian motion, conditions for existence of this 
integral are simpler since assumption ([22|) can be simplified. 
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Theorem 5. Let assumptions (Ai), {A3), {A2), {A'^) {Bi) and {B2) hold and let 
function he non-random and satisfy the following assumption: 

(-B5) There exists such p > Q that 

rp2H-l+p 

lim sup < 00. 

T^oo ip'^{t)dt 

Then estimate is strongly consistent as T ~> 00. 

Proof. It follows from |MMV) and the Holder inequality that for any r > 



E 



Denote Ft 



TT 



An = 



{Fn>£}- 



Jo v>Ht)dt 
Then 



Also, for any > 1 and any e > define event 



E\ (p(s)dB^\'- ll'/'lli 1 [0,JV] N^H-^)r 

P{An) < e-^ . < C{H,r)---^ < C{H,r) 



(/o v^it)dty 



I'^IIl2[o,a'] 



I'^IIl2[0,7V] 



Under condition [Be,) we have P{Am) < C{H,r,p)N 2 . If r > |, then it follows 
immediately from the Borel-Cantelli lemma that series P{An) converges, whence 
Fn — > a.s. as — > 00. Now estimate the residual 



i?jv — sup 

Te[N,N+l] 



Ft — F, 



N\- 



Evidently, 



Rn < sup 



I^^WdB^ 



TelN,N+i] ip'^{t)dt 



and it is sufficient to estimate 



sup 

TelN,N+l] 



J^ip{t)dB^ ^ suprg[jv,Ar+i] 



Fn, 



lN^it)dB," 



:== R 



N- 



According to Theorem 1.10.3 from |MishQ8| and the Holder inequality. 



E 



sup 

Te[Ar,Ar+l] 



T 
N 



^(t)dB^ < C{H,r)M\l ^ < C{H,r)M\l 

' IT 

Now we can use condition (B^ to conclude that for any e > 

II'^IIl2[o,jv] 

< C7(ff,r)e-ni^||Z:[o,^] < C{H^r)e-^ N-^^-"^'^^-^). 

and apply the Borel-Cantelli lemma again. Then we obtain 



We can set r > 



2ff-l+p 



that i?-^ a.s. as A^ 0, which means that 9\. ' is strongly consistent. 

Theorem 6. Let assumptions (Ci) and (C2) hold, and, in addition, 
(C3) ipl{s)ds ^ 00 a.s. 



□ 
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(C4) There exist such a > 1 — H and p > 1 that 
T"+'^-\\ogTrll](Po+^)is)\ds 

Then estimate o!^^ is strongly consistent as T 00. 



(23) Tp — ^ > a.s. as I ~> 00. 



Proof. The last term in the right-hand side of © tends to zero under condition 
(C3). The proof of convergence of the second term repeats the proof of Theorem 

a □ 

Similarly to Theorem[5l conditions stated in Theorem[6]can be simplified in case 
when function is non-random. 

Theorem 7. Let assumptions (Ci) and (C2) hold. Then, if functions ipi and 
are non-random, function ipi satisfies condition (B^), function tp2 is bounded, then 
estimate o!^^ is strongly consistent as T ^ 00. 

Now we shall take a look at the properties of sequential estimates. 

Theorem 8. (a) Let assumptions (-Bi) — (^3) hold. Then estimate ^^^(j,-) is 

unbiased, efficient, strongly consistent, — 0)^ — j^, and for any 

estimate of the form 

with r < 00 a.s. and E JJ" x'id(M-^)s < h we have that 

E{e%~ef <E{er-6f. 

(b) Let function Lp be separated from zero, \^{s)\ > c > a.s. and satisfy the 
assumption: for some 1 — H < a < 1 and p > 

. ^ C'''^^ m+^){s)\ds 

(2) 

as h 00. Then estimate ^^t-(^) is strongly consistent. 

(c) Let function tpi be separated from zero, \^p{s)\ > c > a.s. and let function 
•d satisfy the assumption: for some \ ~ H < a < 1 and p > 

i,^^'^ m+^){s)\ds 

as /i —> 00. Then estimate 'is strongly consistent. 

(d) Let function be non-random, bounded and positive, ipi be separated from 
zero. Then estimate d^^^i^-f is consistent in the following sense: for any 

V 

-5-0 as /i —> 00 . 



p>^,E 



e - 6^'^ 



v{K 
r-r{h) 



Proof, (a) Process /J" ' XsdM^ is a square-integrable martingale which implies 
that estimate is unbiased. Besides, the results from |LipSh| , Chapter 17, can 
be applied to directly, therefore estimate 6*^^^!^^ is efficient, E{9''^^^^ — (?)^ = i. 
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and for any estimate of the form 0^ — °2^^j^j-H^ = + j'x ^a^^jvfg^ with t < oo 

a.s. and E xld{M"), < h we have that £;(6l^^(i) - 0)^ < EiOr - ef. Strong 
consistency is also evident. 

(b) We have that \ J^''^\{s)dB^ \ < {T{h))"+"-^+P J^^'''' \{V^^^){s)\ds. It is 

sufficient to note that h — Jq^^^ fids > c^T{h). The proof of statement (c) is now 
evident. 

(d) It was proved in [MishOSj that in case of non-random bounded positive func- 
tion < i9(s) < "(9*, for any stopping time v 



I ^ t 

E[ sup 



0<t<-u 







d{s)dB^\^^y <C(H,p)r(EvP"y 



Furthermore, same as before, |</?i(s)| > c and h = j^^'^'^^ ip\{s)ds > c^v{h). These 
inequahties together with the Burkholder-Gundy inequahty yield 



E 



9(3) 



□ 



Remark 3. Another proof of statement (a) is contained in [Prara] . Assumptions (p4)) 
and (j25p hold, for example, for bounded and Lipschitz functions Lp and •& correspond- 
ingly. 

4.2. Linear models and strong consistency. I. Consider the linear version of 
model ([5]): 

dXt = 0a{t)Xtdt + b{t)XtdB^, 

where a and b are locally bounded non-random measurable functions. In this case 
solution X exists, is unique and can be presented in the integral form 



Xt = Xq 



6 [ a{s)Xsds+ [ b{s)XsdB^ ^ xoexp{9 [ a{s)ds + [ b{s)dB^\. 
Jo Jq ^ Jo Jo 



Suppose that function b is non-zero and note that in this model 

^^^^ = m- 

Suppose that ip{t) is also locally bounded and consider maximum likelihood estimate 
9!^\ According to (jG)), to guarantee existence of process J', we have to assume that 
the fractional derivative of order j — H for function (;{s) := ip{s)si^^ exists and is 
integrable. The sufficient conditions for the existence of fractional derivatives can 
be found in [SMKj . One of these conditions states: 

(Bq) Functions cp and <r are differentiable and their derivatives are locally inte- 
grable. 

So, the maximum likelihood estimate does not exist for an arbitrary locally bounded 
function if. Suppose that condition (Bq) holds and limit <Jo = linis^.o <r(s) exists. 
In this case, according to Lemma [T] and Remark [2J process J' admits both of the 
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following representations: 



and assuming (B3) also holds true, the estimate 9)p' is strongly consistent. Let us 
formulate some simple conditions sufficient for the strong consistency. The proof is 
obvious and therefore is omitted. 

Lemma 3. If function ip is non-random, locally bounded, satisfies (Bq), limit (^(0) 
exists and one of the following assumptions hold: 

(a) function ip is not identically zero and Lp' is non-negative and non- decreasing; 

(b) derivative <;' preserves the sign and is separated from zero; 

(c) derivative is non- decreasing and has a non-zero limit, 

then the estimate o!^^ is strongly consistent as T 00. 

Example 1. : Let the coefficients are constant, a{s) = a ^ and b(s) = 6 7^ 0, then 
the estimate has a form oip"^ = 6 -\- ^^^^^i^m and is strongly consistent. In this 
case assumption (a) holds. In addition, power functions <^(s) = s^ are appropriate 
ioT p > H — 1: this can be verified directly from ([B]). 

Let us now apply estimate 9)p ' to the same model. It has a form (|5]). We can use 
Theorem [5] directly and under assumption {B^) estimate Oj, is strongly consistent. 
Note that we do not need any assumptions on the smoothness of (p, which is a clear 
advantage of 6*^ . We shall consider two more examples. 

Example 2. : If the coefficients are constant, a(s) — a ^ Q and h{s) = 6 7^ 0, then 

the estimate has a form = 6 + . We can refer to Theorem [5] and conclude 

that 6rp' is strongly consistent. Alternatively, we can use Remark [5] which states 
that \Bi}\ < ^T"{\ogT)P for any p > 1 and some random variable ^, therefore 

(2) (2) 

a.s. as T — ^ 00. In this case both estimates 9)p and 9^ are strongly 
consistent and E{9 — 9]j, ) — ^ ^^(j-i — has the same asymptotic behavior as E[9 — 



Example 3. : If non-random functions and <r are bounded on some fixed interval 
[0, io] but is sufficiently irregular on this interval and has no fractional derivative 
of order | — _ff or higher then we can not even calculate J'{t) on this interval 
and the maximum likelihood estimate does not exist. However, if we assume that 
ip{t) ~ fH-i+p infinity with some p > 0, then assumption (B5) holds and estimate 
9rp' is strongly consistent as T — > 00. In this sense estimate Oj, is more flexible. 

II. Consider a mixed linear model of the form 




(26) 



dXt = Xt{9a{t)dt + b{t)dB^ + c{t)dWt), 



16 



Y. KOZACHENKO*, A. MELNIKOV^* AND Y. MISHURA* 



where a, b and c arc non-random measurable functions. Assume that they are 
locally bounded. In this case solution X for equation (|26p exists, is unique and can 
be presented in the integral form 

Xt=xoexp{e [ a{s)ds+ [ b{s)dB^ + [ c{s)dWs ~ ^ [ c^{s)ds]. 

^ Jo Jo Jo 2 Jo J 

In what follows assume that c(s) ^ 0. We have that (pi{t) — and f2it) ~ 
Estimate Qj. has a form 

^27) g,(3) ^ \l ^MdYs ^ ^ ^ ^x{s)^i{s)dB^ ^ >fr{s)dWs 
lo ^lis)ds ipl{s)ds ipl{s)ds 

In accordance with Theorem [71 assume that function Lpi satisfies [B^) and (^2 is 
bounded. Then estimate Oj, is strongly consistent. Evidently, these assumptions 
hold for the constant coefficients. 

4.3. The fractional Ornstein-Uhlenbeck model and strong consistency. I. 

Consider the fractional Ornstein-Uhlenbeck, or Vasicek, model with non-constant 
coefRcients. It has a form 

dXt = 9{a{t)Xt + b{t))dt + j{t)dB^, t > 0, 

where a, b and 7 are non-random measurable functions. Suppose they are locally 
bounded and 7 = j{t) > 0. The solution for this equation is a Gaussian process 
and has a form 

Xt = e«^W (xo + 9 b{s)e'^^'^'Us + J* ^{s)e''^^'UBf^ E{t) + G{t), 

where A{t) = ^^a{s)ds, E{t) = e^"^(*)(^a:o -I- 9 ^^b{s)e^^^^''^ dsj is a non-random 

function, G{t) — e^"^^*^ ^{s)e~'^ ^^"^ dB^ is a Gaussian process with zero mean. 

Denote c(t) = = Now we shall state the conditions for strong 

consistency of the maximum likelihood estimate. 

Theorem 9. Let functions a, c, d and 7 satisfy the following assumptions: 

{B-j) -ai < a(s) < -02 < 0, -ci < c(s) < -C2 < 0, < 71 < 7(5) < 72, 
functions c and d are continuously differentiate , c' is hounded, c'{s) > 
and c'(s) — >■ as s 00. 

Then estimate 9^^ is strongly consistent as T ~> 00. 

Proof. We shall check the conditions of Proposition [TJ Obviously, il;{t, x) = c{t)x + 
d{t) G Ci(K+) X C2(M) and 

J{t)^ f lH{t,s){d{s) +c{s)E{s))ds+ [ lH{t,s)c{s)G{s)ds := F{t) + H{t). 



Furthermore, assumptions (Ai), (A3), {A'2), (A'^) and (Bi) hold. Note that the 
trajectories of process G are a.s. Holder up to order H, whence 

lim si-"c{s)G{s) = 0. 

Therefore 

J'{t) = F'{t) + H'{t) = F'{t) + / lH{t, s)f{s)G{s)ds + / lH{t, s)c{sh{s)dBf , 

Jo Jo 
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where f{s) — (^^ — H^s^^c{s) + c'{s) + 9a{s)c{s). Evidently, is Gaussian pro- 
cess with mean and variance that are bounded on any bounded interval. There- 
fore, condition {B2) holds. As for condition (-B3), we must verify that Irx, = 
{J[)^t^^~^dt — 00 a.s. For any A > consider the moment generation function 

QtW = Ecxp{-\It} = Ecxp{-X [ {J[ff"-^dt} 

Jo 



and 



/>oo 

600 (A) = ^exp{-A/oo} = ^exp{-A / {J[ff"-^dt}, 

Jo 



so that 9oo(A) — limT^oo 6t(A). Evidently, 



whence 



OriX) < e«(A) := i?exp{ - A( J'^t^-'^dt)"]. 



Random variable Jj-t^^^ dt is Gaussian with mean M{T) and variance a'^{T), 
say. Note that for a Gaussian random variable f = to + <tA^(0, 1) we can easily 
calculate 

(28) E exp{-aC'} = [^aa^ + l) ' cxp { - ^£5^ } • 

This value attains its maximum at the point to = 0. Hence, it is sufficient to prove 
that 

^lim e(?)(A) := ^li^m^i^exp{ - ^{[Kt"-^'} = 0. 

_i 

However, it follows from (^5]) that 0^^(A) = { ^^^"^ ^' ^^^^^^0'"^ prove 

the strong consistency of the maximum likelihood estimate we only need to 
analyze the asymptotic behavior of a^. More specifically, we need to prove that 

^ ^ 00 as r — ?> 00. In what follows we apply the following formulae from |NVV99j 
and |MMV] for Wiener integrals w.r.t. the fractional Brownian motion 

E f ' gis)dB^ f ' h{s)dB^ = H{2H - 1) f f g(si)/i(s2)|si - S2\^"-^dsids2 
Jo Jo Jo Jo 

< CiH)\\g\\L 1 [o,ti]II^IU 1 [o,t2]- 

TT TT 

(a) Let e <0. Divide H[t"-^dt into two parts: J^f H[t"~^dt = H!j}^+H^\ 
where 



and 



H^'> ^ / t"-i I lH{t,s)f{s)G{s)dsdt 



H^^ ^ I t"-^ [ lH{t,s)c{s)j{s)dB^dt. 
Jo Jo 
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Since functions c and 7 are bounded from below and from above, 



Jo Jo Jo Jo 

(29) x\si~S2\^"-^dsids2dtidt2-C{H)f fiht2)"-^ 

Jo JO 

rti rt2 



X 



r r ^^=l,2^H{k,Si)\sl - S2\^"~^dsids2dhdt2 X C {H)T'^ 
JO JO 



as T — >■ 00. 

Consider the behavior of /. Under assumption (Bj) terms s^^c(s) +c'(s) vanish 
at infinity, 6a{s)c{s) is negative and separated from zero. Therefore, there exist 
Ci > 0, i — 1,2 and so > such that — Ci < /(s) < —C2 for all s > so- Bounded- 
ness of / implies that E{H^^)'^ has the same asymptotic behavior as 



{tit2)"-^ r f\u,^i^2lH{U,s,){-f{s,))) 

So J So J So J So 

Si pS2 ^ , fSi 

X 



7(ui)7(u2)exp|6'f 

So -'So ^Jui 

(30) + / ]a{v)dv>\ui ~ U2\^^^^duidu2]dsids2dtidt2 

J U2 

>C{H) f j iht2)"-^ r [\i1,^i,21h{U,s,)) 

J So J So J So J So 

x( f I \ui-U2\^"^^duidu2]dsids2dtidt2>iC{H)T^. 

■J So J So 



Relations ((29)) and (pO| mean that the asymptotic behavior of is x C{H)T^ 
and ^ — ?> 00 as T — 7> 00. 

(b) Let 9 > 0. This case is more involved. The asymptotic behavior of E{H^^)'^ 
is the same as before, C{H)T^, since it does not depend on 6. As for E{H^^)'^, 
denote = lH{t, s)f{s)G{s)ds, then 



(1) 



i.H- 



2dt = T"-^K 



H 



t^-^Kfdt. 
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In addition, denote r{t) = exp{—0 a{s)ds}, tp{t) — f{t) cxp{9 a{s)ds}. Apply- 
ing Pubini theorem several times, we obtain that 



sds 

el. 



= T"-i ^ lH{T,t)f{t)G{t)dt - - ^ ^ lH{t,s)f{s)G{s)dsdt 

= T"-h [ lH{T,t)ii^{t) [ r{s)dB^dt 
Jo Jo 

-(h-^) [ t"-i [ lH{t,u)i){u) [ r{s)dB^dudt 
^ 2/ Jo ^0 Jo 

= r r{s) r lH{T,t)mdtdB^T"-"^- 

Jo Js 

{h-^I^ r{s) t"-i lH{t,u)tp{u)dudtdB^ 



= r{s)(T"-i lH{T,t)i){t)dt 
.{h-^ t"^^ lH{t,u)i^{u)dudt^dB^ . 



Denote 



F{T, s) = ^ Ih{T, t)tp{t)dt -{h-^ J, u)^{u)dudt 

= t-^-"e^ So f(t) {t"--^{T -tf--" 

/ u"-i{u-t)^-"dv^dt:=Fi{T,s)-F2{T,s). 



and 



F+(T,s)=Fi(T,s)+F2(T,s). 



Function / is bounded, positive for s > so and separated from zero. For the sake of 
technical simplicity, we can put f{t) = a{t) = 1. Besides, we can omit the constant 
multiplier cr- Then 



< E{H^^f = / / e^'e'^'F{T,s)F{T,t)\s - t\''"-^dsdt 
Jo Jo 

< / / e'^'e^'F+{T,s)F+{T,t)\s-t\^"-^dsdt. 
Jo Jo 
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Consider the terms containing Fi{T, s)Fi{T,t): 

e^"e^*Fi(T, s)FiiT,t)\s - tf"-^dsdt 

Jo Jo Js 

X j^^ vi-"{T - v)i-"e-^^dv\s - t\^"-^dsdt 

Jo Jo Js 

X j\T - v)^-"e-'^^^-*^dv\s - t\^"-^dsdt. 

Applying Holder inequality we conclude that integral fJ{T — u)i~^ e~^^"~^^ du 
admits the following bound: 

J s 

<2"-iiT-s)i~" ' e-'^^-'^du+[j^^^{T-uf-^"duY[j^^^e-'^^-^Uuy 
< C{H) ((T - s)3-« + (T - sf-^y 



Therefore 



h<C{H)(T^"-^ J J {st)i-"{{T-t){T-s))i-"\s-t\''"-^dsdt 



+Tj J {st)^-"\s-t\'''-'dsdtj <C{H)T\ 
Furthermore, function e^^F2{T, s) admits the following bounds: 

e'^'F2{t,s)<C{H)si-" t"-i{T - t)i-"e-^^*-'Ut 

J s 

J s 

Note that function jj t^~ie~^^*~^^dt decreases in s since its derivative equals 
e^%jjt"-ie-^'dt -s"-i)<0. Therefore, 

e^'F2{t,s) < C{H)Ti-"si-" C t"-'^e-'^'dt < C{H)Ti-" s^-" . 

Jo 
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The latter implies that the term containing F2{T, s)F2{T,t) admits the following 
bounds: 



l2 = r r e''e'*F2{T,s)F2{T,t)\s-t\^"-^dsdt 
Jo Jo 

<C{H)T^-^" / {st)i-"\s-t\'^"-^dsdt<CiH)T 
Jo Jo 



So, E{H^^^ f X C{H)T*''^" asymptotically and if we compare this to asymptotical 

behavior of E{h!^^Y x C{H)T^, we can conclude that ^ x C{H)T^ ^ oo as 
T ^oo. 

(c) Let 6 = 0. Then it is easy to verify that E{H^^ f x C{H)T and we can 
refer to the case 9 > {). □ 

Remark 4. The assumptions of the theorem are fulfilled, for example, if a(s) = —1, 
h{s) = h € M. and 7(3) = 7 > 0. In this case we deal with a standard Ornstein- 
Uhlcnbeck process X with constant coefficients that satisfies the equation 

dXt = 9{b - Xt)dt + 7dBf , t > 0. 

This model with constant coefficients was studied in [KlLeBr] where the Laplace 
transform 9t(A) was calculated explicitly and strong consistency of o!^^ was estab- 
lished. Therefore, our results generalize the statement of strong consistency to the 
case of variable coefficients. 

IL Consider a simple version of the Ornstein-Uhlenbeck model where a = 7 = 1, 
6 = xq = 0. The SDE has a form dXt = OXtdt + dB^^ , t > with evident solution 
Xt = e^* e~^^dB^ . Let us construct an estimate which is a modification of : 



^e-O'dB^ 



2 



9(2) ^ ^ ^ 



~(2) 

Theorem 10. Let 9 > 0. Then estimate 9rp is strongly consistent as T 00. 
Proof. Applying Remark [S] yields 

\j\-''dBf\ < e-'^lB^l+J^-'^lB^lds < ^{e-''^T"+^+j\-'\s"+Pds) < C, 

where is a random variable independent of T. So, it is sufficient to establish that 
/o°° ( lo ds = to prove the strong consistency of 9!^K Similarly to the 

proof of Theorem [9l we can consider the moment generation function 

Eexp{-\J (y e-^"dBf)^ds} <£;exp{-Ar-i(y" J e-'^'^dSf ds) ^} 
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where 



([ r e-^^dB^dsY ^ f I r f e-'^-'^\u-v\^"-^dudvdsdt 
^ Jo Jo ' Jo Jo Jo Jo 

Jo Jo Jo Jo 

2 



Jo Jo Jo 
1 



^r^2HQ-2^j (^l_e-«-^jdsj xr2«0-2, 

whence the proof follows. □ 

Appendix A. 

To apply Theorem [2] to the fractional derivative of the fractional Brownian mo- 
tion and to prove Theorem[3l we need an auxiliary result. In what follows we denote 
by C{H,a) a constant depending only on H and a and not on other parameters. 

Lemma 4. Let Zi > for i = 1,2. In addition, letO<H<l,l — H<a<l and 

1^ ^ \2H .,2H .2H 
I = Z22(^+"-l) + 2^2(if+a-l) ^ \Z2 - Zl| - Z, - Z, 



Then I < C{H,a)\z2 - 

Proof. Let Z2 > zi > (the case zi > Z2 > can be dealt with in a similar way). 
We can rewrite / as 

/ = (zf +"-1 - zf +"-1)2 + 2(ziZ2)^+"-l 

+ ((^2 - Zlf" - {Z^ - Zff - 2{z,Z2f ){z^Z2r-^ 

_ r^H+a-l H+a-l-.2 , (z2 - Zi)^ - (z2 - Zi Y _ j , t 
-{Z2 -Z, ) + ^— ^^-^^ -/1+/2. 

Recall a simple inequality }f — < (6 — aY for 6 > a, < r < 1. Since < 
iJ + a — 1 < 1, we can estimate Ii by (z2 — zi)^^-^^"-^^ . Furthermore, I2 can be 
rewritten as 

(ZiZ2)^ "(Z2 - Zi)^l.-"+" 

where « = f > 1, /(«) = ^-X-m-^-^^ ^ 

Calculate the limit of function / at 1: 

(u-l)2g-(^^-l)2 
hm /(.) = hm (^„^^2(^+.^i) ■ 



Here 

lir,, 



and 

hm , ^"'!27ffy D = 1™ - 1)'-'^^'" = 0, 
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since lim - — ^ — H. Calculate the limit of the function f at infinity: 

< lim = lim Jfe = 0. 

This implies that function / is bounded, i.e. there exists C{H^ a) > such that 

l2<C{H,a)iz2-zif^''+"-'\ 
and the proof follows if we combine the bounds for 7i and l2- □ 

We are now ready to check conditions (-D2) and (D3) for the fractional derivative 
of the fractional Brownian motion. 



Lemma 5. Let 

X{t) = ^ + / ^du, 

where < t2 < ti, Q < H < 1,1- H < a < I. 

Then the following hounds hold: 

1 ) for any < ^2 < ii 

{E{X{t)fy^ < C{H,a){ti-t2f 

2) (a) Let H + a < |. Then for any < t2 < ti, < S2 < si and any 
0<e<(ff + a-l)Ai 

{E\Xit)-Xis)\^)i 

< C{H,a){l + - sil V |t2 - S2|))^+"-'-^(ii V s,))' 

with C{H,a) not depending on X, its arguments and e. 

(b) Let H + a> ^. Then for any < t2 <ti, 0< S2 < si 

{E\X{t) - X{s)f )i < C{H,a){\ti - si| V \t2 - S2\)Hti V si)"+"-i. 



Proof. The first statement follows immediately from the Minkowski's integral in- 
equality: 



In order to prove the second statement, denote X\{t) = \ and X2{t) 



III Evidently, 

(31) (£|X(t) - X(s)|2)^ < {E\X^(S) - Xi(s)|2)5 + {E\X2{t) - X2(s)|2)i 
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Let ti > si, the opposite case can be considered in a similar way. Then 

(£;|Xi(t)-Xi(s)p)^ 

/ / BP - BP 

= E 



{ti-t2y-" (<i-s2)i-" ih-s^y-" (si-s2)i-° 



(32) 



< E 



Bll Bii - B,^ 



(tl-i2)l-" (tl-S2)l- 



/ / - Bf -Bi , , 

It is more convenient to estimate the squares (/s)^ and (14)^ from ([5^ instead of 
I3 and I4. As for (/s)^, we can calculate it explicitly and then estimate it with the 
help of Lemma m {I^y can be evaluated similarly. 

EiB^-B»yB"-Bg) 



(tl-<2)l-"(il-S2)l-" 

2 



(il-<2)l-"(tl-S2)l-" 



(33) X [tr - ^ (^r + - ih - t^y") - \ {tr + 4" (ti - S2y") 

+\ {ti" + 4" - 1<2 - .^2^)] = {h - 12)^(^+"-^) + (ii - s,y<^"+--^) 



\t2 - 82?" ~ {h - t2y" - (ti ~ S2y" , _ ,2iH+a-l) 

(tl-t2)l-"(il-S2)l-" <C-l^,"j|t2 S2\ 



We derive from (1551) that 

(34) /3<C(i/,a)|i2-S2|''+"~\ 
and similarly, 

(35) /4 < C(i7, a) |ti -si|^+"-\ 
It follows immediately from ([34]) and ([SS]) that 

(36) (£;|Xi(t) - Xys)\^)-2 < C{H,a) (|ti - sil V ^2 - 521)^+""' 
Now estimate 



t2 Js2 (^-S2)2-"'^"j I 



(/ rti tdH tdH nsi r>H tdI 
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Let, for instance, < ^2 < S2 < si < ii (other types of relation between these 
points can be handled similarly). Then 



S2 nH _ nH \ 2\ 2 



t2 (u-h) 



2\ 2 



Using the Minkowski's integral inequality we immediately obtain 



(38) 



2\ 2 



H + a — 1 



Similarly, 



(39) ^^^i7T7^(^-^^)""""- 



Again, using the Minkowski's integral inequality and Lemma |4] we conclude that 
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(40) 



h < 



E 



- Bg B^ - Bl 



+ 



(w-i2)^-" (U-S2)2-" 
{S2-t2?'' -iu-hf" 



du 



(u-i2)2-«(w-S2)2-« 



du 



+ 



{S2-t2f" -{u-t2?" -{U-S2?" 



{u-t2Y-°'{u- S2) 



1-a 



du 



J S9. 



,2H 



du 



'du 



<C{H,a) f \u-S2)-Hu-t2)-Hs2-t2)"+"-h 

J S2 

+C{H,a) f \u-S2)"+"-\u-t2)-Hs2-t2)Uu=:Is+l9. 

J S2 

Evidently, 

/8 = {S2 - t2)"+"-' f\u - S2)---{U - t2)---du = {S2 - t2f+"-'ll0 

up to the constant multiplier and for any < £ < ^ integral Iio can be rewritten 



as 



'^du 



ho= {U- S2)~^{u-t2)~^ 

J S2 

= / {y + 1) ^dy< [- 
Jo V- 



< 



Sl - S2 



S2 - t2 

1 1 



+ ^y 

{y + l)-^y-^-'dy < C{l + e-^) ('^-'^ 



,52-^2 

Therefore, for any 0<£<(iJ + Q!-l)Ai 

(41) h < C{H, a)(l + £-i)(s2 - t2)"+"-^-%si - S2r. 



S2 - t2 
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Furthermore, 



h = (S2 - ^2)5 / {U - S2f+'^-\u - t2)—-du = (S2 - 12)^ hi 



up to a constant multiplier. In the case when H -\- a < ^ the integral In can be 
rewritten as 



/ll= r {U~ S2)"+''-^{u-t2)-"^du 

J S2 



' S2 



JQ 

In case when H + a > ^ integral In admits an obvious bound 

hi < r (u - S2)"+"-^iu - 82)-^ du < C{H, - S2)^+"-' . 

Finally, ioi H -\- a = | integral In admits the same bound as Iiq. Therefore, 

(42) l9<C{H,a){s2-t2)"+''-' 
for + a < |, 

(43) /g < CiH,a){s2-t2)Hsi~S2)"+"-^ 
for H + a > I , and 

(44) l9<CiH,a){s2~t2)^-%si~S2Y 

for iJ + a = I . 
This implies that 

(45) ^^(t,s) < C{H,a){l + e-'){\ti - si\ V 1^2 - S2\)"+''-'-' {si V hY 

for iJ + a < |. In case iJ + a > | we can put e = + a - | e (0, 5) in (|41]) and 
conclude that 

(46) ^^(t,s) < C{H,a){\ti - si\ V \t2 - S2|)^(si V ti)^+"~i. 

The proof follows immediately from ((3T|) and p6p - (j46|) . □ 



Proof of Theorem[3l First of all we should verify conditions (Di) — (-D3). Condi- 
tion (Di) is evident, since X is continuous in both variables. According to the 2nd 
statement of Theorem[51 condition {D2) holds with /3 — e, < e < {H + a — 1) A ^ 
and 7 = _ff + a — 1 — ein case when a + H < | , and with (3 = H + a — ^ and 7=5 
in case when a + H > | . According to the first statement of Theorem [SJ condition 
(D3) holds with S = H + a-l. 

Let A{t) = {t"+°'-^\\ogt\P) V 1 for some p > 1 and for any i > and let 
bi^e\l>0. Then 61 = {e"-"+°'-^hP) V 1 and A{bi) = e'(^+"-i'. Therefore, in 
this case series S{S) converges since 

p(i+l)(H+a-l) °° 

SiS) = + J: = e«--(l + m < 00. 

i=i 1=1 
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Moreover, it is easy to check that 1 + ^ — ;^ = for any values of a + if, hence 
Ki = 0. This imphes that all conditions of Theorem [5] hold true and we can apply 
the theorem with A{t) = logij^") V 1 which concludes the proof. □ 

Remark 5. Instead of the fractional derivative, we can consider the fractional Brow- 
nian motion itself and apply the same reasoning to it. This case is much simpler 
and we immediately obtain that supQ<g<;( < {{t^ {log{t)y) V l)C(j') for any 
p> 1. 
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